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Abstract 

The  research  performed  imder  the  above  titled  ONR  grant  has  resulted  in  a  number  of 
important  advances  in  the  treatment  of  integrodifferential  treinsport  equations,  such  as  those 
obtained  when  the  Fokker-PIanck  coefficients  depend  on  the  “temperatmre”  of  the  particle 
distribution,  as  defined  by  a  suitable  integral.  In  sximmary,  we  have  developed  a  new  technique 
for  independently  determining  the  temperature  as  a  function  of  time  using  the  energy  moments 
of  the  initizd  distribution.  This  is  equivalent  to  decoupling  the  original  integrodifferential  problem 
into  two  parts:  (i)  determination  of  the  temperature  variation,  and  (ii)  determination  of  the 
particle  distribution  as  a  function  of  time  and  energy.  This  represents  a  major  advance 
in  the  treatment  of  integrodifferential  Fokker-PIanck  equations  because  the  two 
problems  obtedned  can  be  solved  separately  without  recourse  to  traditional  predictor- 
corrector  alogorithms.  The  technique  has  been  applied  with  success  to  the  problem  of  self- 
consistent  Comptonization  in  astrophysical  plasmas.  The  theoretical  approeich  utilizes  a  novel 
combination  of  physics,  applied  mathematics,  and  advanced  symbolic  computational  techniques, 
and  shows  good  promise  for  generalization  to  the  treatment  of  other  integrodifferential  systems  of 
interest  to  the  Navy.  The  multi-disciplinary  CSI  program  at  George  Mason  University  has  proven 
an  ideal  venue  for  this  investigation. 
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I.  Introduction 

During  reidiation-dominated  X-ray  transients  in  AGNs,  energy  is  rapidly  transferred 
between  the  photons  and  the  electrons  until  the  electron  temperatinre  (Te)  has  equilibrated  to  the 
inverse-Compton  temperature  of  the  radiation  (Tic)-  This  occurs  on  a  much  shorter  timescale 
than  the  equilibration  of  the  radiation  into  a  Wien  distribution.  Once  the  temperatures  have 
equilibrated,  the  condition  Tg  =  Tic  continues  to  be  enforced  by  the  Compton  exchange  of  energy 
between  the  photons  and  the  electrons.  The  radiative  tremsfer  problem  then  becomes  highly 
nonlinear,  since  the  variation  of  the  electron  temperature  determines  the  spectrmn  (through  the 
Kompaneets  equation),  while  the  spectnun  determines  the  electron  temperature  (through  the 
condition  Te  =  Tic).  It  is  therefore  extremely  useful  to  develop  techniques  for  independently 
calculating  the  electron  temperature  as  a  function  of  time  for  an  arbitrary  initial  radiation 
spectrum. 

In  the  reseztrch  carried  forth  here,  we  have  developed  a  general  technique  for  determining 
the  electron  temperatmre  by  using  the  energy  moments  of  the  initial  spectnun  to  construct  a 
continued  fraction  with  dramatic  convergence  properties.  The  electron  temperatmre  so  obtained 
is  then  folded  back  into  the  PDE  to  compute  the  spectral  solution.  Hence  the  original  integrodifferential 
problem  is  decoupled  into  two  parts:  (i)  determination  of  the  electron  temperature  as  a  function 
of  time,  and  (ii)  determination  of  the  photon  spectrum  as  a  function  of  time  and  energy. 


II.  Time-Dependent  Comptonization 


In  rawiiation-dominated,  fully-ionized  plasmas,  photon  creation  and  destruction  are 
unable  to  establish  locaJ  thermodynamic  equilibrium,  and  the  photons  and  electrons  interact 
primarily  via  Compton  scattering.  Neglecting  thermal  opacity  sources,  the  temporal  evolution 
of  a  distribution  of  photons  in  am  infinite,  homogeneous  plasma  with  electron  number  density  Ug 
and  electron  temperature  Tg  is  governed  by  the  Kompaneets  (1957)  Fokker-Planck  equation 


dn  _  rieO'^h  1  d  ’  4  /_  kTe  dn\ 

dt  mgC  1/2  du  ^  \  h  dv )  ' 


(1) 


where  n[y,  t)  is  the  photon  occupation  number,  1/  is  the  photon  frequency,  eind  the  terms  on 
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the  right-hand  side  express  the  effects  of  .recoil  losses  and  energy  diffusion,  respectively.  Note 
that  the  properties  of  the  electrons  need  not  be  constant,  and  as  we  argue  below,  the  electron 
temperature  will  track  the  inverse-Compton  temperatmre  of  the  photons  during  radiation- 
dominated  flzires. 

A  Maxwellian  distribution  of  electrons  will  exchange  no  net  energy  with  the  photons  if 


Te  =  Tic,  where 

_  1  J^hu^ndv 


(2) 


is  the  inverse-Compton  temperature  of  the  radiation,  defined  here  without  regard  to  stimulated 
processes.  In  a  radiation-dominated  plasma,  energy  is  rapidly  transferred  between  the  photons 
2uid  the  electrons  until  Te  =  Tic-  This  occurs  on  a  much  shorter  timescale  than  the  equilibration 
of  the  radiation  into  a  Wien  distribution  n  oc  exp{—hu/kTe).  If  the  cooling  timescale  is  much 
shorter  than  the  diffusion  timescale  (td  »  tc)>  then  most  of  the  spectral  evolution  takes  place 
after  the  electron  and  inverse-Compton  temperatures  have  equilibrated.  Dmring  this  later  phase, 
the  condition  Te  =  Tic  is  enforced  by  the  Compton  exchange  of  energy  between  the  two  species 
(Becker  &  Begelman  1986),  and  Compton  scattering  drives  the  photon  distribution  toward  a 
Wien  distribution. 

The  physical  requirement  that  Te  =  Tjc  in  the  Kompaneets  equation  increases 
the  complexity  of  the  problem,  since  the  equation  becomes  integrodiflferentieJ  in 
nature.  The  temperature  must  therefore  be  computed  self-consistently  along  with  the  spectrum, 
starting  from  a  specified  initial  radiation  distribution.  The  coupled  integrodifferential  system 
is  usually  solved  numerically  using  predictor-corrector  techniques.  However,  we  show  here 
that  the  variation  of  the  electron  temperature  can  be  determined  in  advance  using 
only  the  initial  spectrum  itself.  With  the  electron  temperature  independently  determined, 
the  solution  for  the  spectrum  can  be  generated  straightforwardly  by  numerical  integration  of  the 
Kompaneets  equation.  In  §  III  we  develop  a  new  technique  for  computing  the  perturbation  series 
for  the  temperature  using  only  the  energy  moments  of  the  initial  spectnun.  In  §  IV  the  resulting 
(divergent)  power  series  is  transformed  into  a  (convergent)  continued  fraction  and  evaJuated  for  a 
variety  of  initial  spectra. 
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III.  Perturbation  Method 

The  Kompaneets  equation  can  be  rewritten  in  terms  of  the  dimensionless  energy 


TUeCr 

and  the  Compton  y-parameter  (essentially  the  dimensionless  time) 

/•* 


JO 


dn  _  \  ^  \  ^  ^  ^ 

dy  dx  L  W  * 


where 


._kTeit) 
m,<?  ' 

Next  we  introduce  the  energy  moments  of  the  radiation  spectrum,  defined  as 


h{y)=  f 
Jo 


ndx , 


and  note  that  the  photon  energy  and  number  densities  can  be  expressed  as 


Upl,  =  ^  J  hv^ndv  =  Sirme<^  (‘^)  ‘ 


We  also  note  that  the  condition  Tg  =  r,c  requires  that 

4/3’ 

and  operate  on  equation  (5)  with  x^  dx  to  obtain  the  differential  recurrence  relation 


{l+l)h- 


4/3-f«+l 


Since  there  can  be  no  net  energy  exchange  between  the  photons  and  the  electrons  (Te  =  T,,,), 
it  follows  that  h  =constant.  Equations  (10)  and  (11)  can  therefore  be  applied  recursively  to 
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obtadn  the  y-derivatives  of  9  in  terms  of  the  moments  of  the  radiation  spectnnn.  In  particular, 
the  derivatives  of  9  evaluated  at  y  =  0  can  be  expressed  in  terms  of  the  moments  of  the  initial 


radiation  distribution,  which  are  known.  For  example,  the  zeroth  and  first 
y  =  0  are  respectively 


derivatives  of  9  at 


(12) 

(13) 


etc.  In  Figure  1  we  plot  the  Taylor  series  for  9(y)  obtained  in  the  case  of  a  bremsstrahlung  initial 
spectrum, 


n(x,y)Uo  =  x-3e-"/"®°. 


(14) 


The  series  is  obviously  divergent  beyond  a  small  radius  of  convergence.  However,  we  show  in  §  IV 
that  the  series  can  be  recast  as  a  convergent  continued  firaction. 

IV.  Continued-Fraction  Representation 


Although  the  perturbation  series  obtained  using  the  energy  moments  of  the  initial 
spectrum  is  often  qmte  divergent,  useful  information  about  the  temperature  variation  can 
nonetheless  be  recovered  by  using  the  power-series  coefficients  to  construct  the  continued-fraction 
representation 


^{y)  = 


1  1 

c^y 

i  1 

n-- 

c^y 

(15) 


l-l-_ 

l  +  CpfV 

which  is  simply  a  rationed  function  approximation  to  the  unknown  function  9{y).  The  continued- 
fraction  coefficients  cq,...,c^  are  computed  using  the  algorithm  derived  by  Becker  (1988).  The 
resulting  sequence  of  continued  fractions  is  plotted  in  Figmre  2  for  the  case  of  a  bremsstrahlung 
initial  spectnun.  In  contrast  to  the  highly  divergent  Taylor  series,  the  sequence  of  continued 
fractions  is  uniformly  convergent  to  a  unique  limit  function.  The  general  form  of  the  temperature 
variation  is  in  good  agreement  with  the  numerical  calculations  and  semi-analytic  models  of 
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Becker  &  Begelman  (1986).  Furthermore,  the  nm  of  d{y)  computed  by  substituting  the  resulting 
solution  for  the  photon  spectrum  n{x,y)  into  equations  (7)  and  (10)  agrees  nearly  perfectly  with 
the  limit  function  depicted  in  Figme  2. 

We  have  also  obtained  results  for  monochromatic  initial  spectra  and  for  combinations  of 
power-laws  and  exponentials,  and  we  find  that  the  sequence  of  firactions  is  convergent  in  all  cases, 
although  imiform  convergence  is  achieved  only  when  the  initieil  photon  number  spectrum  diverges 
logarithmically  at  low  energy.  The  observed  X-ray  spectrum  for  fast  (unresolved)  transients  can 
be  computed  by  averaging  n{x,y)  over  the  characteristic  diffusion  time  for  the  photons. 

V.  Summary 

•  We  have  developed  a  new  technique  for  independently  determining  the  self-consistent 
variation  of  the  electron  temperature  in  Comptonizing  plasmas. 

•  The  method  assumes  that  Te  =  at  all  times,  and  utilizes  only  the  information 
contained  in  the  energy  moments  of  the  initial  spectrum. 

•  The  sequence  of  continued  fractions  obtained  in  the  case  of  a  bremsstrahlimg  initial 
spectrum  demonstrates  uniform  convergence  to  a  unique  limit  function  9{y),  which  is  the  correct 
temperature  variation. 

•  The  technique  essentially  decouples  the  integral  and  diflferential  parts  of  the 
original  integrodiflferential  Kompaneets  equation,  thereby  eliminating  the  need  to  use 
a  conventional  predictor-corrector  algorithm. 

•  The  method  cem  be  generalized  to  treat  other  integrodiflferential  equations  of 


special  interest  to  the  Navy. 
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Fig.  1.  —  Sequence  of  Taylor  series  approximations  of 
order  0  through  24  for  the  temperature  function  theta [y] 
for  a  bremsstrahlung  initial  spectrrm. 
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Fig.  2.  -•>  Sequence  of  continued-fraction  approximations 
of  order  0  through  24  for  the  temperature  function  theta [y] 
for  a  bremsstrahlung  initial  spectrum. 


